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Abstract 

The statistical description of oscillator strengths for systems like hydro- 
gen in a magnetic field is developed by using the supermatrix nonlinear 
cr-model. The correlator of oscillator strengths is found to have a univer- 
sal parametric and frequency dependence, and its analytical expression 
is given. This universal expression applies to quantum chaotic systems 

with the same generality as Wigner-Dyson statistics. 
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A conventional way to approach the phenomenon of quantum chaos is to study 
spectra of quantum systems These spectra are known to be well described by 
Wigner-Dyson statistics of energy levels 0, based on the eigenvalue distribution of 
random matrix ensembles. Usually chaos manifests itself in highly excited states of a 
system. Excitation from a regular state (e.g., the ground state) to a chaotic state can 
be achieved by applying an external perturbation. Experiments of this type provide 
independently the transition probabilities from a given initial state to excited states, 
as well as the spectrum of the system. These probabilities for optical excitations 
are known as oscillator strengths. Like the eigenenergies, the oscillator strengths of 
highly excited states are random, and should be described statistically. Along this 
line, the distribution of the transition matrix elements was predicted to have the 
Porter-Thomas distribution in the chaotic limit j3|. 

In this Letter we develop a statistical description of oscillator strengths based 
on the Wigner-Dyson hypothesis. We investigate a statistical characteristic different 
from Ref. || — the correlation function of oscillator strengths at different frequencies 
and "magnetic fields". It is emphasized that we make no further assumptions such 
as on the form of the matrix elements. Since the information contained within the 
oscillator strengths is independent of the spectral properties, their statistics provide 
a characterization of quantum chaos which complements the usual energy level corre- 
lations. We determine universal parametric correlations of oscillator strengths which 
should apply with the same generality as Wigner-Dyson statistics. We have also dis- 
covered an unexpected differential relation between this parametric correlations and 
the parametric density of states correlator of Ref. j3|. 

The results of this Letter can be applied to responses of various impurity states 
in semiconductors, or to microwave and optical excitations of chaotic quantum dots 
etc. They can be easily modified to describe a much broader class of phenomena, e.g., 
cross sections of nuclear excitations by neutrons. However, as a concrete example, we 
have in mind the hydrogen atom in a magnetic field, which attracts great interest as a 
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"real" quantum system with a chaotic classical analogue. The energy spectrum of this 
system is known to obey Wigner-Dyson statistics ||. In addition, the distribution of 
the transition matrix elements is observed to obey the Porter-Thomas distribution ||. 
The oscillator strength per unit energy is conventionally defined as 

W(Q, B) = 2053 1($ | z \a)\ 2 5{E a {B) - E Q {B) - ft), (1) 

a 

where B denotes the magnetic field. More generally, B can stand for any external 
parameter of the system. In writing down Eq. ([!]), we use the transition matrix 
element between a certain initial state |$o) with the energy Eq and a highly excited 
state I a) with E a induced by the light polarized along z-axis. However, it is not 
necessary to assume this particular form of the transition matrix elements to proceed 
our theory. 

We have evaluated the correlation function of oscillator strengths 

s(n,B) = (w(n,B)w(n + n,B + B)), (2) 

for small energy separation (ft <C ft) and at slightly different "magnetic fields" (B -C 
B). Here (• • •) denotes the averaging over a certain interval of the spectrum and/or 
of the parameter B. The magnitude of the correlator S(Q, B) is not universal. It 
depends on the initial state |$o) ; on the type of excitation, and on the direction of 
the polarization. (See Eq. ( |TTD below.) However, the dependence of S on ft and B 
turns out to be determined by a universal function s(co,x), 

S(n,B) = S [s(u,x) + l}; S = (W) 2 , (3) 

where we use the rescaled parameters defined by M 

u= ft/A ; A = (E a+1 - E a ) , (4) 

X=B^/C{0); C(0) = ((dE a /dB) 2 )/A 2 . (5) 

The two terms in Eq. (^) correspond to the connected and disconnected parts of 
the correlator. The function s(u, x) is classified according to the Dyson universality 
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classes: orthogonal (o), unitary (u), or symplectic (s). It is convenient to present 
s(lj, x) by 



( d 2 d \ 

s e (uj, x) = I — + 2Ce^ J Ma;,x), 



(6) 



where e = o, u, s; ( a = ( s = 2, ( u = 1 and 
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In Eq. (H), the integral region is defined by A G (—1,1), Ai )2 G (1, oo) for the orthog- 
onal case, or A G (1, oo), Ai G (—1, 1), and A 2 G (0, 1) for the symplectic case. 

Let us compare S(Q, B) with the more conventional statistical characteristic of 
quantum chaotic systems. For a certain value of B, the density of states at energy E 
can be written as v(E, B) = J2 a ~ E a (B)). The parametric correlator of densities 
of states is written in terms of rescaled variables Eqs. (|4],[|) 

(9) 



u(E, B) v(E + tt, B + B)) = — [1 + k(u, x)] 



The explicit forms for k D (uj, 0), k u (u, 0) and k s (u, 0) were obtained by Dyson (see, e.g., 
Ref. ||), and were recently evaluated at finite x [Q by the supermatrix method 0. 
In our notation these results can be written as Hi 



k(u,x) = d 2 h(uj, x)/du 2 . 



(10) 



Both k(u,x) and s(u,x) are singular at u = and x = 0. We present in Fig. 1 the 
ratio r(u>,x) = (1 + s(w,x))/(l + k(u,x)) which is not singular. 

Comparing s(u, x) of Eq. with k(u), x) of Eq. QT0[) , we get a differential relation: 



duj- 



;S(U,X) 



2C 



_d_ 

dx" 2 



k(ui, x). 
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This relation is not trivial at all || , especially when considered from the semiclassical 
point of view [|1],|9[]. That approach is based on the approximation of the density of 
states v(E, B) by a certain sum over the classical periodic trajectories. Although 
the averaged oscillator strength itself may be expressed by the summation of the 
periodic trajectories [|10||, the oscillator strength correlator apparently requires all 



closed trajectories to be taken into account, not just the classical periodic ones. (See 
also Ref. [JTTJ] for some discussions.) Eqs. (|6|JlOD for s(uj,x) and k(uj,x) are valid for 



arbitrary uo and x, including nonperturbative regime at to ?S 1. We believe that these 
formulas as well as Eq. (|TTD have wider region of applicability than the conventional 
semiclassical approach. 

It is worth mentioning that the definition of h(u, x) should be interpreted rather 
formally, and in fact, integral expressions Eq. (0,§) for h(u),x) diverge logarithmically. 
Of course, the corresponding integrals for s(u, x) and k(u, x) are converged. However, 
we have introduced h(uj,x) not only to make formulas for s(u),x) and k(u,x) more 
compact. Note that the finite part of h(uj,x) — h(0,x) is related to the generalized 
parametric number variance function v(u,x) by (compare with Ref. \T2\) 



h(uj, x) - h(0, x) = - (v(u, x) - v(0, x)) , (12) 

where v(uj, x) = [n(e + u, x + x) — n(e, x) — oj} 2 and n(e, x) = J2i #( e ~ e i( x ))- 

Now let us sketch the derivation of our main results (fj-||). One can write S(fl, B), 

5(0, B) = Jdr 1 dr 2 dr 3 dr 4 A(r 1 , r 2 , r 3 , r 4 )i2(ri, r 2 , r 3 , r 4 ), (13) 



where 
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A = —z^zsz^Kn^o^mrs)^^, (14) 

7T Z 



R = (imG* + ^( ri , r 2 ; B) ImG| o+n+n (r 3 , r 4 ; B + B)) . (15) 



and A (r, r'; B) = (r\{E - H(B) ± iO)" 1 ^') 

are the retarded and advanced Green's 
functions. The factor A is determined by the initial state wave function $o( r ), whereas 



R reflects only the properties of excited states, which are presumed chaotic. In 
contrast to the evaluation of K(Q,B) in Eq. (|9]), we now have to accommodate the 
dependence on four different positions ri, r 2 , r 3 ,r 4 as depicted in Fig. 2a. 
We apply the zero-dimensional supermatrix model |7| to determine R: 

^R(rx, r 2 , r 3 , r 4 ) = /(r 12 )/ (r 34 ) Re (l - QllQH 



-/(r u )/(r M ) Re (QIIQH) q - /(r 14 )/(r 32 ) Re (QHQH) q , (16) 
where the dimensionless function f(r) = (lmG R (r)) / (lmG R (0)) which is not uni- 



versal |14j]. For example, /(r) becomes sin(pr) exp(— r/2£)/pr in three dimen- 
sional metallic grain with the Fermi momentum p and the mean free path £. 
(■ • -)n of Eq. (|T6| ) denotes the integration over the saddle point manifold (• • -)n = 
JDQ (• • •) exp (-(?7rcj/4)STrQA - (7r 2 x 2 /32C e )STr[Q, A] 2 ), and the definition of STr 
as well as A and the structure of supermatrices Q can be found in [0. We can 
identify physically each term of Eq. ( |16"D as in Fig. 2(b-d), - - DOS fluctuation- 
type (Fig. 2b), Cooperon-type (Fig. 2c) and Diffuson-type (Fig. 2d). In fact, we 
get the DOS correlator of Eq. (H) by integrating Fig. 2b over ri and r 3 , whereas 
integrating Fig. 2c and 2d over ri and r 3 (Fig. 2c) or r l5 r 2 (Fig. 2d) gives another 
function which was recently shown to correspond to the dielectric function of peri- 
odic systems with chaotic unit cell fl3| . Definite integration over Q-matrix leads to 
Re(QllQH - l) Q = -2 - 2d 2 h(u,x)/duj 2 , Re(QllQH} Q = -4 dh(uj , x) / dx 2 , and the 
Cooperon-type contribution {Q1 2 QH)q is equal to the Diffuson-type (QHQIDq for the 
orthogonal and symplectic cases, or vanishes for the unitary. Everything combined, 
we obtain our results (§-§), with the prefactor S equal to 

-l 2 



So 



2Qu / c?r 1 cir 2 2; 1 2 2 <l>*(r 1 )/(r 12 )<l>o(r 2 ) 



(17) 



since $o( r ) is real for orthogonal and symplectic cases. As it was mentioned, the 
constant Sq as well as the parameters A and C(0) is not universal. However, the 
dependence of S on the rescaled frequency uj and "magnetic field" x is universal. 

6 



In conclusion, we have developed the statistical description of oscillator strengths 
of a chaotic system such as hydrogen in a magnetic field. In particular, we have ob- 
tained the analytical form of the correlator of oscillator strengths S(Q, B), dependent 
on Q and B (Eqs. ([|fj)). We have found an identity Eq. ( |TTD that relates S with the 
density of states correlator. It looks like a challenge for the semiclassical theory of 
quantum chaos to derive this differential relation. The universality of S(Q, B) is as 
general as that of the correlation function of densities of states K(Q, B), and should 
be applicable to various physical situations: impurity states in semiconductors, to 
microwave and optical excitations in chaotic quantum dots, etc. We believe that ex- 
perimental studies of the statistics of the oscillator strength provides a new way of 
characterizing the chaotic systems, complimentary to the spectral statistics. 
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FIGURES 



Orthogonal 




FIG. 1. Schematic plot of 1 + s(cj, x) normalize by 1 + k(co, x) for orthogonal (above) 
and unitary (below) cases. 
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